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ABSTRACT 

We present a wide-field imaging study of the colors of bright galaxies (< M*+2) in a sample of 
12 X-ray selected clusters and groups of galaxies at z ~ 0.3. The systems cover one of the largest 
ranges in X-ray luminosity (Lx ~ 10 43 - 10 45 erg s _1 ), and hence mass, of any sample studied at this 
redsliift. We find that the ‘red’ galaxies form a tight color-magnitude relation (CMR) and that neither 
the slope nor zero-point of this relation changes significantly over the factor of 100 in X-ray luminosity 
covered by our sample. Using stellar population synthesis models we find our data allow a maximum 
possible change of 2 Gyrs in the typical age of the ‘red’ galaxies on the CMR over the range of Lx of 
our sample. We also measure the fraction of blue galaxies (f&) relative to the CMR in our clusters and 
find a low value of fj, ~ 0.1 consistent with other X-ray selected cluster samples. We find that there is 
no correlation between f;, and Lx over our large Lx range. However, we do find that both the CMR 
and ft, do depend significantly on cluster radius, with the zero-point of the CMR shifting blueward in 
B — R by 0.10 ± 0.036 magnitudes out to a radius of 0.75 times the virial radius. This color change is 
equivalent to a luminosity weighted age gradient of ~ 2.5 Gyrs per log(radius) and is consistent with 
previous studies of the radial change in the zero-point of the CMR. It thus appears that the global 
cluster environment, in the form of cluster mass (L_y), has little influence on the properties of the 
bright cluster galaxies, whereas the local environment, in the form of galaxy density (radius), has a 
strong effect. The range of ~ 100 in Lx corresponds to a factor of ~ 40 in ram-pressure efficiency, 
thus suggesting that ram-pressure stripping, or other mechanisms that depend on cluster mass like 
tidal stripping or harassment, are unlikely to be solely responsible for changing the galaxy population 
from the ‘blue’ star forming galaxies, that dominate low density environments, to the ‘red’ passive 
galaxies that dominate cluster cores. 

Subject headings: galaxies: clusters: general - galaxies: elliptical and lenticular, cD - galaxies: 
evolution — galaxies: fundamental parameters 


1. INTRODUCTION 

The first clear evidence for galaxy evolution in clus¬ 
ters was presented by Butcher and Oemler (BO) in a 
series of papers in the late 70s and early 80s. Initially, 
they demonstrated that the fraction of blue galaxies 
(f h) in two high redshift c entrally concentrated clusters 
(IButcher fc Oemlei1ll978al) was signi ficantly higher than 
that found in similar local clusters ( IButcher V Oemler! 
fi978bl) . This suggests that unlike the cores of lo¬ 
cal compact clusters which contain almost exclusively 
quiescent early-type galaxies, similar clusters at high 
redshift contain many star forming galaxies in addi¬ 
tion to the early-types. This result was subject to 
significant scrutiny, particularly the required statis¬ 
tical removal of foreground and background ga l axies 

S ieu fe Spinradl 119811 iDressler fc Gunnlll982t iKod 
In an attempt to resolve this, IButcher fc Oemlerl 
studied 33 galaxy clusters up to a redshift of 
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0.5. This showed an increasing trend of blue frac¬ 
tion with redshift, but with a large scatter. Several 
othe r studies have reported similar trends. For instance 
[Rakos fc Schombert! ( 1995 1 find an even faster increase, 
but do not limit the ir clusters to cen t rally con centrate d 
ones, and Margoniner. de Carvalho. Gal Kr. Piorgovskil 
( 200 11 also confirm a similar trend using a large sam¬ 
ple of Abell clusters. 

Follow-up observations of the BO blue galaxies 
showed spectroscopic signatures of galaxies with ei- 
ther ongoing or recen t ly-truncated star form ation 
(iCouch fc Sharpie^ 119871 IDressler fc Gunnl I1983H , and 
morphological studies with HST have shown most 
of these galaxies to be late-type spirals but with 
an increased incid ence of a distur bed or interact- 
ing appear ance (IDre ssler. Oem er. Butcher. & Gunn) 

S; 

ICouch et alJ 


Couch. Ellis. Sharole^faSmaiT 


199 


However, i t is not clear t h at the increased b l ue fra ction 
is universal. iSmail. Edge. Ellis, fc Blandfordl (119981) . us¬ 
ing 10 bright X-ray clusters at z ~0.25, have found no sig- 
nif icant ev i dence fo r an increase in the blu e frac tion, and 
lAnington-Smith. Ellis. Zirbel. Kr. Oemlerl (ll 9931) showed 
no increase in radio selected groups. It has also been sug¬ 
gested that the clusters studied at high redshift are not 
similar to those th a t were observed locally. For example, 
lAndreon Kr. Ettorii ( 1999 ) have shown that the Butcher 
and Oemler clusters have increasing X-ray luminosities 
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with redshift. 

Whilst evidence exists for strong evolution in the 
late-type population of clusters, the early-type galaxies 
have shown very little evidence of evolution beyond the 
passive ageing of their stellar populations. In particular, 
studies of the correlation betw een the co l ors and mag- 


1978a H) in clusters (|Bower. 

Laicev. & Ellis! 1992albl: 

Terlevich. Caldwell. & Bower 

200l[) have shown 


that this relation has a small sc atter and e volves 
very little to high redshift dKllis et all 1997 1 
[Stanford. Eisenhardt. fc Dickinsonl Il998t1 . This slow 
evolution in the slope, zero-point and scatter of the 
color-magnitude relation (CMR) suggests that the stellar- 
populations of early-type galaxies were formed at high 
redshift (z > 2), evolving purely passively since, with 
the slope being consistent with a metallicity seq uence 
rather than an age s equence llKodama fc Arimotoill997t: 
[Kodama et al.lll99il. 


The only evolution in the early-type population in clus¬ 
ters appears to be in the number of SO galaxies. A spec¬ 
troscopic and HST study (MORPHS) of 10 rich galaxy 
clusters (0.37 < z < 0.56) found that the fraction of ellip¬ 
tical galaxies remained constant, b ut that the f raction of 
SO galaxies decreased with redshift llDressler et al.ll997ll . 
Further, they found an increase in both the star forma¬ 
tion rate and in the num ber of galaxies show ing evidence 
of recent star formation iPoggianti et al.lfl999t1 with red¬ 
shift. However, a similar spectroscopic study of high red¬ 
shift rich clusters, this time X-ray selected, found a very 
weak redshift depende nce of the sta r formation rate with 
redshift ( e.g. CNOC: IBalogh et al.llll999f) . 

The decrease in SO galaxies and the increase in blue 
star forming galaxies and galaxies with recently trun¬ 
cated star formation with redshift have often been linked. 
It has been suggested that spiral galaxies falling into 
the cluster experience either a burst of star formation 
followed by its rapid truncation or just a rapid trunca¬ 
tion of their star formation. This results in the fading 
of the disk and the formation of an SO galaxy. Several 
mechanisms that can occur in the cluster environment 
may be driving this transformation. These include 
ram-pressure str ipping EE nji A Gottlll972l IFuiitalll998t 
Fjyitk^^JaggghhngL 1999; Abadi. Moore, fe Bower 

19991 lFuiitall2001t iQuilis. Moorev&^owe r 2000), galaxy 

h arassment (|Moore_e^dJ 19961 [Moore. Lake, fc Kata 
119981 IMoore. Lake Quinn, fc Stade]|ll999H t idal eff ects 


(Bvrd fc Valtonen 199(3; iHenriksen & Bvrd 1996| Fuiita 

1998 

), 

and strangulation llLarson. Tinslev. & Caldwel 

1980 

2001 

IBalogh. Navarro, fc Morris 2000; iDiaferio et al 

)• 

These processes have the affect of removing an 


infalling galaxies gas supply and thus halting ongoing 
star formation. This idea has been further reinforced 
by the change in galaxy properties with cluster ra dius . 
The Morpholo gy- Radius/Density relation l|Oemlerl|1 9741 
iDressleil 1 1980 1 at low redshift is an excellent example 
of this and first inspired the idea that a galaxy’s star 
formation may be halted by the cluster e nvironment re¬ 
sulting in a morpho logical transformation llGunn fc Gotti 
ir97llOemleHri97l . T his relation was further shown to 
exist at high redshift llDressler et all 119971 iTreu et all 
2003). 

Similarly, [Butcher fc Oemler ( |1984l> showed that the 
blue fraction increased with radius in their clusters and 


this tren d has generally been confirmed by other stud¬ 
ies (e.g. Ab rahgunet al.: 1.191)fi IKod ama fcZBpwd 12QQ1; 
iFairlev et alJ l2002f ). iBalogh et all <1199711 show that 

the star formation rate depends strongly on cluster ra¬ 
dius in the CNOC sample. In the same sample, they 
find no radial dependence of the fraction of post-star- 
burst /star-forming galaxies , unlike the MORPHS clus¬ 
ters (IPoggianti et 23 EH). Both, however, conclude 
that the truncation of star formation is independent of 
the morphology-density relation. Recently, two large 
studies using the latest galaxy red shift surveys (2dF, 
ll.ewis et a 1 .1121)021 iGdmez et a.l .1120031 SDSS) have found 
very clear relationships between star formation and den¬ 
sity in the vicinity of clusters. Both studies find a char¬ 
acteristic density (radius) at which the average star for¬ 
mation rate rapidly drops. 

Some observational evidence also exists for an en¬ 
vironmental depende nce of the prope rties of elliptical 
galaxies in clusters. IA bra,h a,m et a.l J i ll 996f i find that 
the galaxies on the CMR become bluer as a func- 
tion of cluster r adius in A2390 at a redshift of 0.23. 
[Terlevich. Caldwell, fc Boweij j|2001]( find a si milar trend 
in the Coma cluster, as does lPimhblet, et' in a 

combined sample of 11 rich clusters at z~0.1. 

When trying to draw general conclusions regarding the 
evolution of cluster galaxies as discussed above, one par¬ 
ticular concern is how cluster selection may have affected 
the results. It has already been noted that the BO sam¬ 
ple has systematic differences as a function of redshift. A 
further concern comes from the fact that the clusters used 
were those that already had available dat a. This meant 
that they were often unusual objects llSandage et alJ 
EH, which was the reason they had been studied 
originally. This ty pe of selection has often occ u rred 
(e.q. MORPHS; (Stanford. Eisenhardt. fc Dickinsonl 
ttaas) and when other more careful selection criteria 
are used, differing results ar e fou nd {e.g. CNOC; 
iSnrail. Edge. Ellis, fc Blandfordl 1199,S i. It should also 
be noted that the clusters studied have mainly been 
massive rich clusters, as they provide the best statistics 
and are easiest to locate at high redshift. This may 
be of particular concern as they are very rare objects 
and are unlikely to be the evolutionary predecessors 
of the clust ers with which they are compared locally 
( Kauffmann 1995 1. Galaxy properties clearly depend 
on environment, and it may be that by studying only 
the most massive clusters a non-representative picture 
of galaxy evolution has been generated. 

This paper presents further results from a comprehen¬ 
sive study of galaxy populations in X-ray selected clus¬ 
ters at both intermediate and high redshift. The clusters 
were selected to cover a wide range in X-ray luminosity 
(Lx) and all were imaged over a wide field in three bands. 
The wide range in X-ray luminosity and hence mass and 
the wide-field imaging allows us to simultaneously in¬ 
vestigate how the cluster galaxy population d epends on 
both cluster mass and cluster radius. In IFairlev et alJ 
iE?M Paper I) we presented measurements of the blue 
fraction in an initial sample of 8 clusters over a large 
redshift range. In this paper 12 intermediate redshift 
clusters from the sample are used to investigate how the 
blue fraction and CMR depend on both X-ray luminosity 
and cluster radius. 

Throughout this paper a cosmology with Dm = 0.3, 
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Oa = 0.7 and Hq = 70 km s 1 Mpc 1 is assumed. 

2. CLUSTER SELECTION AND DATA 

The 12 clusters presented here are selected to cover 
a large range of X-ray luminosities and hence cluster 
masses over a small redshift interval at z ~ 0.3. Because 
they are all selected using their X-ray emission and are at 
similar redshifts, comparison of their galaxy populations 
should be free of any selection effects driven by the clus¬ 
ter selection technique. The clusters were drawn from 
4 X-ray selected cluster sa mples: t he Exten ded Medium 
Sensitivity Survey fEMSS: IOioia et a!~ll 1 99(1) . the South¬ 
ern Serendipito us High-redshift Ar c hival ROSAT Clus - 
ter (SSHARC; iCollins et alJ Il99?t iBurke et al.l l200.il) . 
the Bright Serendi pitous High-redshift Archival ROSAT 
Cluster ('BSHARC lNichol et alll999tlRomer et al .1200111 
and the Wide Angle ROSAT Pointed S urvey (WARPS; 
iScharf et alJ 119971 iPerlman et a. 1.1 1200211 . The low flux 
limits of SHARC and WARPS allow the selection of clus¬ 
ters with low Lx, whereas the large area of the EMSS 
allows the selection of rarer high Lx clusters. Table ^ 
lists the name, original catalog, coordinates of the cen¬ 
ter, redshift, X-ray luminosity, and telescope (with which 
the observations were made) of the cluster sample. 

2.1. Observations 

Each cluster was imaged over a wide field of view 
through three filters, B , V and R, in the Johnson- 
Cousins system. These filters correspond approximately 
to U, B , and V at rest at the cluster redshift, minimiz¬ 
ing any K-corrections required. The imaging was un¬ 
dertaken using two telescopes: the Issac Newton Tele¬ 
scope (INT) at the Roque de Los Muchaclios Observa¬ 
tory on the Island of La Palma and the Blanco Tele¬ 
scope at the Cerro Tololo Inter-American Observatory 
(CTIO), Chile. Each telescope was equipped with a 
wide-held mosaic CCD camera. Th e IN T Wide Field 
Camera lllves. Tulloch. Churchillll996t) consists of 4 
thinned EEV 2kx4k CCDs mounted at the INT prime 
focus in an ‘L’ configuration. Each CCD has a pixel size 
of 13.5 /cm corresponding to 0.33 arcseconds on the sky. 
There are 2048 by 4100 pixels in the useful imaging area 
of each CCD, giving a field of view of 22.8 by 11.4 ar- 
cminutes. The 4 CCDs together give a total imaging area 
of ~ 0.29 deg 2 . The CTIO MOSAIC II camera consists 
of 8 2048 by 4096 SITe CCDs mounted in 2 rows of four. 
Each pixel is 15 /im corresponding to 0.27 arcseconds on 
the sky, giving each CCD dimensions of 9.2 by 18.4 ar- 
cminutes, and a total imaging area of ~ 0.38 deg 2 . The 
gaps between the CCDs are ~ 40 pixels or 11 arcseconds. 

The INT data were acquired on 4 nights from 17-20th 
June 1999. The observing conditions were generally good 
through the 4 nights, except for the first night during 
which light cirrus cloud was present. Despite this night 
being non-photometric, observations of clusters were still 
made using slightly longer exposure times with calibra¬ 
tion images of the same fields taken on the following 
night. It was also noted that Saharan dust was present 
in the upper atmosphere throughout the observations. 
This is a common feature at the Roque de Los Mucha- 
chos Observatory and is known to show ‘gray’ absorp¬ 
tion, increasing the atmosp heric extinction equally over 
all wavelengths llKinell 198511 . The CTIO data were taken 
on the nights of the 26th and 27th of September 2000. 


The observing conditions were good on both nights with 
no cloud cover. Tables [2] and 0 give details of the obser¬ 
vations. 

2.2. Data Reduction 

The basic reduction ste ps were undertaken using the 
IRAF package MSCRED ( Va.lded Pl998lf . This package 
was written specifically to reduce data from the two 
NOAO mosaic cameras, one of which is the MOSAIC 
II camera on the CTIO Blanco telescope. The main ad¬ 
vantage of this over other data reduction packages is its 
ability to operate on the multi-CCD data simultaneously, 
without the need to split each image frame into sepa¬ 
rate CCD images and process them individually. It was 
also possible to use this package with non-NOAO Mosaic 
data such as from the INT WFC. Before this could occur, 
modifications to the headers of the INT images describ¬ 
ing the relative position of the CCDs to each other were 
required. 

Both data sets were reduced in a similar manner. Only 
when a different procedure was applied to each data set 
will the individual data sets be identified. 

In general the data reduction proceeded in the stan¬ 
dard manner for CCD imaging. After bias subtraction 
the images were flat-fielded using dome flats. The R band 
data from the INT had to be corrected for fringing. A 
fringe frame was generated from the object frames and 
then subtracted from the individual frames. Sky flats 
were then generated by median combining the object 
frames, and a further flat-field correction applied using 
these. Typically this produced frames with a maximum 
variation in the sky level of < 0.5% and < 1% for the 
CTIO and INT images respectively. Cosmic rays, sat¬ 
urated pixels and stellar bleed trails were identified and 
added to the bad pixels masks. In addition to these steps 
a linearity correction was required for three of the INT 
WFC CCDs. This non-linearity was investigated previ¬ 
ously by the Cambridge Astro nomic al Survey Wide Field 
Survey team iMcMahon et al.ll999tl and found to be con¬ 
stant with time. They characterized the non-linearity 
and have provided equations to convert the measured 
values of CCDs 2, 3 and 4 to the value of CCD 1. A fur¬ 
ther correction was also required for the CTIO MOSAIC 
camera data to remove the ’cross-talk’ in the controller 
electronics. This can cause pixel values for one ampli¬ 
fier to be affected by another and induce ‘ghost images’. 
This is well understood for the CTIO MOSAIC camera 
and a model has been developed to correct the data. 

Wide-held imaging requires a step not normally nec¬ 
essary in CCD data reduction. The curvature of the 
focal plane causes a significant variation to the pixel size 
across the CCD array. After flat-fielding all pixels have 
the same sky level which is incorrect for pixels of differ¬ 
ent sizes and leads to a varying photometric zero-point. 
To correct for this an accurate astrometric solution was 
determined for each image, which was then used to rebin 
the image to a constant pixel scale. Before rebining, the 
pixel scale changed by ~ 4% in the INT WFC and ~ 8% 
in the CTIO MOSAIC II camera from the center to the 
edge of the array and would thus cause a significant error 
in the photometry across the field. 

Before combining the individual images into final 
stacked images of each cluster in each band, the pho¬ 
tometric scales were matched, correcting for variations 
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TABLE 1 

Details of the cluster sample 


Cluster 

Catalogue RA Dec 

Redshift 

Lx 

Lx error 

Telescope 

RX1633.6+5714 

WARPS 

248.4254 57.2366 

0.239 

0.10 

0.02 

INT 

RX0333.0-3914 

SSHARC 

53.2755 -39.2305 

0.245 

0.12 

0.02 

CTIO 

R.X0210.4-3929 

WARPS 

32.6112 -39.4908 

0.273 

0.16 

0.02 

CTIO 

R.X0054.0-2824 

WARPS 

13.5150 -28.4039 

0.293 

0.24 

0.03 

CTIO 

RX1606.7+2329 

WARPS 

241.6825 23.4875 

0.310 

0.69 

0.06 

INT 

RX2237.0-1516 

BSHARC 

339.2504 -15.2766 

0.299 

1.09 

0.08 

CTIO 

MS0407.2-7123 

EMSS 

61.6908 -71.2750 

0.229 

1.87 

0.41 

CTIO 

RX1418.5+2510 

WARPS 

214.6279 25.1813 

0.294 

2.96 

0.20 

INT 

RX0256.5+0006 

BSHARC 

44.1370 0.1032 

0.360 

3.44 

0.17 

CTIO 

MS0353.6-3642 

EMSS 

58.8733 -36.5641 

0.320 

7.05 

1.38 

CTIO 

MS1455.0+2232 

EMSS 

224.3129 22.3427 

0.259 

12.81 

0.64 

INT 

MS2137.3-2353 

EMSS 

325.0633 -23.6608 

0.313 

15.48 

2.11 

CTIO 


Details of the cluster sample, catalogue, coordinates of the center, redshift, X-ray luminosity, error in the X-ray luminosity and Telescope. 
Lx is given for the ROSAT band (0.5 - 2.0 keV) in units of 10 44 erg s _1 . 


TABLE 2 

Description of the INT observations. 


Cluster Filter Night Exposure Time(s) Number Seeing 


MS1455.0 

B 

2 

2250 

3 

1.22 


V 

1 

2400 

2 

1.58 


R 

1 

2700 

3 

1.56 

RX1418.5 

B 

3 

3600 

3 

1.23 


V 

3 

1200 

1 

0.93 


R 

3 

1200 

1 

0.87 

RX1633.6 

B 

3 

1200 

1 

1.15 


B 

4 

1800 

2 

1.20 


V 

4 

1200 

1 

0.88 


R 

4 

1200 

1 

0.79 

RX1606.7 

B 

4 

3600 

3 

1.11 


V 

4 

1200 

1 

1.11 


R 

4 

1200 

1 

0.93 


Summary of the cluster observations made at the INT from 17 - 
‘20 th June 1999. For each cluster the filter, the night, the total 
exposure time, the number of exposures and the mean seeing are 
listed. 

in the sky level, extinction and exposure times. This was 
done by comparing the flux of a number of stars that 
appear on each image. The corrected images are then 
stacked using a mean, excluding pixels in the bad pixel 
mask. 

The photometric calib ration was ach ieved using images 
of standard star fields (lLandoltjlll992 ) taken throughout 
each night. Fields covering a large area were chosen 
so that it was possible to place standard stars on ev¬ 
ery CCD with one pointing. For the CTIO data there 
were sufficient standards to fit for a zero-point, an extinc¬ 
tion and a color correction term on each CCD on each 
night. As would be expected the extinction terms were 
all equivalent within the errors and so a weighted mean 
was taken. The zero-point and color terms were then re¬ 
fitted. Although the color terms may be expected to vary 
slightly between CCDs the small size of these corrections 
meant that they were equivalent within their errors and 
a weighted mean was again used to give a final color term 
for each CCD. The zero-points were then refitted using 
the mean extinction and color terms. Again they were 
found to be sufficiently consistent (rms < 0.01) to com¬ 
bine and use as one term for the whole mosaic. Although 
using one photometric solution for the whole array may 
slightly increase the photometric errors it greatly simpli¬ 
fies the photometric measurements for the stacked im- 


TABLE 3 

Description of the CTIO observations. 


Cluster Filter Night Exposure Time(s) Number Seeing 


RX0054.0 

B 

1 

1650 

3 

1.23 


V 

1 

900 

3 

1.12 


R 

1 

630 

3 

1.17 

RX0210.4 

B 

1 

1650 

3 

1.26 


V 

1 

900 

3 

1.25 


R 

1 

630 

3 

1.26 

RX0256.5 

B 

2 

1650 

3 

1.70 


V 

1 

900 

3 

0.87 


R 

1 

630 

3 

0.84 

RX0333.0 

B 

1 

1650 

3 

0.99 


V 

1 

900 

3 

0.86 


R 

1 

630 

3 

0.92 

MS0353.6 

B 

2 

1650 

3 

1.67 


V 

2 

900 

3 

1.49 


R 

2 

630 

3 

1.49 

MS0407.2 

B 

2 

1500 

3 

1.95 


V 

2 

750 

3 

2.06 


R 

2 

630 

3 

1.84 

MS2137.3 

B 

1 

1650 

3 

1.41 


V 

1 

900 

3 

1.21 


R 

1 

630 

3 

1.13 

RX2237.0 

B 

1 

1650 

3 

1.20 


V 

1 

900 

3 

1.07 


R 

1 

630 

3 

1.05 


Summary of the cluster observations made at the CTIO from 26 - 

27 th September 2000. For each cluster the filter, the night, the total 
exposure time, the number of exposures and the mean seeing are 
listed. 

ages since spatial offsets between the individual images 
mean that a single patch of the stacked image can have 
contributions from several CCDs. 

The INT calibration proceed in a similar manner. 
However, fewer images of standard fields were observed 
over too low a range in airmass preventing an accurate 
determination of the extinction. Instead, extinction coef¬ 
ficients were used from the Carlsberg Meridian Telescope 
a photometric monitoring telescope on the observatory 
site. As for the CTIO, data zero-points and color terms 
were first determined for each individual CCD to ensure 
they were consistent and then combined to determine 
zero-points for the whole array. 

2.3. Constructing Galaxy Catalogues 












Environmental dependence of galaxy colors 


5 



16 18 20 22 18 20 22 

Total V magnitude Total R magnitude 

Fig. 1.— The concentration parameter as a function of magni¬ 
tude in the V and R filters. Top panels show cluster MS2137.3 and 
the bottom show RX0256.5. The lines show the division selected 
to divide stars from galaxies. 

Object detection and photometry were performed 
using the Source Ex tractor package (SExtractor; 
iBertin fc ArnoutJl 19961) . Objects were detected on the 
R band images only, with photometry then being per¬ 
formed on all three images using the objects from the R 
detection. A Gaussian filter of size equal to the image 
seeing was used during the detection process to better 
differentiate between faint sources and noise peaks. The 
areas around very bright sources were also masked as 
these can effect SExtractor’s background determination 
leading to poor photometry. B oth total p hotometry, us¬ 
ing a modified Kron aperture ( Kro 3GIII, and aperture 
photometry were performed. 

It is known that SExtractor’s total photometry esti¬ 
mate suffers from a decrease in the included flux with 
decreasing total flux. In order to calculate a correction to 
the total magnitude, isolated galaxies were identified on 
each image. Using the IRAF task PHOT, the magnitudes 
of these isolated galaxies within a large aperture were 
calculated in order to estimate the true total magnitude. 
For bright objects, the difference between the Sextractor 
and true total magnitudes was approximately constant 
at ~ 0.04 magnitudes or ~ 4%. This rose to ~ 0.2 mag¬ 
nitudes or ~ 15% for the faintest objects. These correc¬ 
tions were applied to all the total magnitudes. Apertures 
with a physical diameter of 20 Kpc were used to calculate 
the colors, representing a minimum size of at least twice 
the seeing, thus ensuring that the small seeing variations 
between images do not effect the colors. The observed 
colors of each object were used to apply the color cor¬ 
rections determined from the standard solutions giving 
the final magnitudes for each object. In order to remove 
any spurious detections, objects that were not detected 
in at least two bands were removed from the catalogue. 
Finally, the number of bad pixels present in an object 
was calculated from the exposure masks and added to 
the catalogue. Any object with more than 10% bad pix¬ 
els in a band had the magnitude in that band flagged. 
The final catalogues of all measured objects in a field 
contained between 15433 and 39977 objects. 

2.4. Star-Galaxy Separation 


In order to differentiate stars from galax- 
ies a met h od s imil ar to t hat presented in 
iMetcalfe. Shanks. Fong, fe Jonedl (I1991T) was used. 
A concentration parameter was defined as the magni¬ 
tude within 4 pixels divided by the total magnitude. 
Figure ^ shows the concentration parameter as a 
function of magnitude for two example cluster fields. 
The stars are visible as the horizontal linear feature in 
the plots and are clearly separated from the galaxies 
at bright magnitudes. The solid lines define the region 
chosen to represent the stars in each plot. These criteria 
had to be selected for each cluster individually as this 
parameter is dependant on the image seeing, with 
increasing seeing moving the locus of stars to higher 
values and greater width. The seeing was observed to 
vary across the wide-held of our images, up to 10% in 
the CTIO images and 20% in the INT images. The effect 
of this radial variation in the seeing is to just broaden 
the stellar locus, reducing the magnitude with which the 
separation can occur, but not causing any systematic 
change in the separation with radius. The B-band data 
were not used for the final separation as they displayed 
a larger scatter in the stellar locus. The separation 
was typically effective to an R magnitude of ~ 22. 
Fainter than this, the number counts are dominated by 
galaxies. Any remaining stars will be included in the 
held region in equal proportions as in the cluster regions 
and thus will be removed by the statistical background 
subtraction (see Section l3~2ll . 

To conhrm the reliability of this method, the star and 
galaxy classihcations in the held of RX0256.5 were com¬ 
pared to those o f the Sloan Digital Sky Survey (SDSS; 
lYork et alfeOOfi) which also images this held. The SDSS 
star-galaxy separation is reliable to m/j ~ 22, compara¬ 
ble to the magnitude at which the separation was pos¬ 
sible here. Of the 2558 galaxies, 52 were classihed as 
stars by the SDSS, and of the 1155 stars, 87 were clas¬ 
sihed as galaxies by SDSS. The SDSS star-galaxy sep¬ 
aration is app r oxim ately 95% accurate at these depths 
llLuoton et al.ll200Hi , suggesting that the separation here 
is at least as accurate. 

The hnal step in the construction of the galaxy 
catalogues is to correct for galactic extinction. 
To apply the re quire d corr ectio ns the maps of 
iSchlegel. Finkbeiner. fc Davis! ('1 9981) were used. Typi¬ 
cally, one extinction correction is applied per held. How¬ 
ever, the resolution of these maps is sufficient to map 
variations over the wide helds of the images used here. 
Investigation showed that the extinction could vary by 
up to 25% across a held. Therefore, the extinction was 
calculated for each galaxy individually by using the in- 
terp olation options avail a ble in the code s provided with 
the ISchlegel. Finkbeiner. fc Davis! (I1998H maps. 

2.5. Photometric Checks and Completeness 

Before using the photometry for our scientihc investi¬ 
gation it is reassuring to check its accuracy against ex¬ 
ternal sources and to determine if the observations have 
achieved the desired depth and accuracy. Fortunately, 
two of the helds imaged here have existing photometry 
from other sources that may be used to check the relia- 
b ility of the photome try presented here. 

iTvson et al.l : 1 998 i measured Johnson B and B — R 
photometry for 99 stellar objects in the held of RX1418.5. 
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Fig. 2. — The R band number d ensity-magnitude relation for all 12 cl uster fields. The solid line show s the fitted relation of 
Metcalfe, Shanks, Fong, & Jones (1991). The dashed line shows the relation of Smail, Hogg, Yan, & Cohen (1995). 


Sixty-five of these objects have photometry measured by 
SExtractor, of which 41 are stars. Comparing with our 
B and R photometry for these 41 stars shows a median 
difference of less than 0.03 and a rms < 0.1 showing a 
good consistency between the data sets to m^ = 22 . 

The field of RX0256.5 has photometry available from 
the Sloan Digital Sky Survey in the Sloan filters 
( - u, q, r , i.z). Using the conversions from (Smith et all 
6 2002 ) and the Sloan PSF magnitudes, the equivalent 
total B, V , and R magnitudes were calculated for all the 
stars present in both data sets. The conversions were 
found to be less reliable at either end of the g — r colour 
range and so the sample was restricted to 0 . 5<<7 — r< 
1.2 leaving 1076 stars in common. The median offset 
was found to be < 0.03 in all of the filters with an rms 
< 0.09 to m# = 22, showing excellent agreement between 
the two data sets. 

Figure [3 shows the number magnitude relations 
(NMR) in the R-band using the SExtractor total mag¬ 
nitudes. Reassuringly, the relations are similar for all 
the cluster fields but do show some scatter as would be 
expected from cosmic variance. The mean density can 
vary by almost a factor of two between cluster fields. 
This illustrates the importance of individual estimates 
of the background density when correcting cluster galaxy 
counts (see Section n~21l . P lotted o_n t he R -band NMR are 
t he re lations derived bylMetcalfe^_Shanks. Fon g, fc Jonesl 
(I1991H and ISmail. Hogg. Yan. & Cohenl (Il995n ~ These 
are consistent with the NMRs of the cluster fields pro¬ 
viding another confirmation of the photometry. 

These relations can be used to estimate the complete¬ 
ness of the galaxy catalogues. The NMRs clearly begin 
to turn over as the completeness declines. Table 0 fists 
the estimated completeness in the detection bands and 
the original target equivalent to My = -19 for an ellipti¬ 


cal galaxy at redshift zero. These values were estimated 
by observing the position at which the NMR appears to 
turn over. This level of precision was adequate, as the 
questions considered in this paper typically use galaxies 
at least ~ 1 magnitude brighter than these limits. 

3. ANALYSIS AND RESULTS 
3.1. Defining a Characteristic Radius 

The variation in the mass of our cluster sample means 
that each cluster has a different physical size. In or¬ 
der to compare the galaxy populations between clusters 
we must define some characteristic radius that is equiv¬ 
alent for all the clusters. If we just use a fixed physi¬ 
cal scale it may include just the core of the most mas¬ 
sive clusters and substantia lly beyond the cor e of the 
least massive. To this end [Butcher fc Oemler ( 11978b ) 
defined a characteristic radius, R 30 , within which they 
measured the blue fraction. R 30 is defined as the ra¬ 
dius containing 30% of the total number of galaxies in 
a cluster. To calculate this, it is necessary to integrate 
the cluster density profile up to the edge of the clus¬ 
ter (R 100 ) after subtracting the background galaxy den¬ 
sity. Most other studies of the Butcher-Oemler effect 
have attempted to calculate Ron for their clusters (e.q. 

ISmail Edge. Ellis, fc Blandfordlll998HKodama fc Bowed 

120011 iFairlev et alJl2002t IPimbblet et alJ 120021) . Apart 

from the original Butcher and Oemler papers, details of 
the calculation of R 30 are scarce, particularly of how R 100 
and the background level were determined. 

Here we attempted to determine R 30 by the following 
procedure. For each cluster, the density profile was cal¬ 
culated in 0.25 areminute bins, by counting all galaxies 
whose centers lay in each bin up to a limiting R magni¬ 
tude of 22.5. The area in each radial bin was then calcu¬ 
lated using the exposure masks of each image. Only the 























Environmental dependence of galaxy colors 


7 


TABLE 4 

Cluster X-ray properties. 


Cluster Redshift Lx ROSAT fcT Ry Ry 


RX1633.6 

0.239 

RX0333.0 

0.245 

RX0210.4 

0.273 

RX0054.0 

0.293 

RX1606.7 

0.310 

RX2237.0 

0.299 

MS0407.2 

0.229 

RX1418.5 

0.294 

RX0256.5 

0.360 

MS0353.6 

0.320 

MS1455.0 

0.259 

MS2137.3 

0.313 


10 44 erg s 1 10 44 erg s 1 keV Mpc arcmins 

0.10 ± 0.02 0.24 ± 0.05 1.53 ± 0.27 1.17 ± 0.10 5.16 ± 0.44 

0.12 ± 0.02 0.30 ± 0.06 1.67 ± 0.27 1.22 ± 0.11 5.27 ± 0.43 

0.16 ± 0.02 0.41 ± 0.06 1.86 ± 0.28 1.26 ± 0.10 5.02 ± 0.40 

0.24 ± 0.03 0.66 ± 0.10 2.19 ± 0.31 1.34 ± 0.10 5.10 ± 0.38 

0.69 ± 0.06 2.24 ± 0.24 3.47 ± 0.33 1.66 ± 0.08 6.08 ± 0.29 

1.09 ± 0.08 3.85 ± 0.35 4.28 ± 0.31 1.86 ± 0.07 6.98 ± 0.26 

1.87 ± 0.41 3.63 ± 0.95 4.32 ± 0.52 1.98 ± 0.12 9.03 ± 0.54 

2.96 ± 0.20 12.47 ± 1.10 6.69 ± 0.45 2.34 ± 0.08 8.87 ± 0.30 

3.44 ± 0.17 14.71 ± 1.01 6.93 ± 0.48 2.25 ± 0.08 7.46 ± 0.27 

7.05 ± 1.38 17.53 ± 4.31 8.13 ± 0.36 2.52 ± 0.06 9.03 ± 0.22 

12.81 ± 0.40 36.74 ± 1.69 5.60 ± 0.29 2.20 ± 0.06 9.15 ± 0.25 

15.48 ± 2.11 45.65 ± 8.42 5.00 ± 1.10 1.99 ± 0.22 7.23 ± 0.80 


The calculated bolometric luminosities, temperatures and virial radii for all clusters i] sing the T.y-T relatio n (>f [VikhlhiijLj^tjdJ |200^). 
For MS0353.6, MS1455.0 and MS2137.3 the temperatures quoted are those taken from IMushol.zkv Sell aril i 1 99A I and Labia n 

ITOPOL The errors in Ry are calculated from the Lx error along with the errors on the Lx-T relation. 




0 5 10 15 


Radius / arcmin 


Fig. 3.— The total number of galaxies as a function of radius 
from each cluster’s center, for a representative sub-sample of our 
clusters. Typical error bars are shown above the lines. 

areas which had maximum exposure in all three bands 
and the galaxies within those areas were included. A 
further correction to the area was required to take into 
account the effects of large objects and crowding. Large 
objects potentially obscure galaxies in front of or behind 
them. This can be a problem, particularly in the cores of 
dense clusters where many bright ellipticals lie. For all 
objects larger than 900 pixels, an estimate was made of 
the area brighter than the deblending threshold. These 
estimates were then summed and subtracted from the 
previously calculated area. 

Figure [3] shows the integrated number of galaxies as 
a function of radius for four example clusters. The in¬ 
tegrated number was calculated by subtracting a back¬ 
ground density typically calculated as the mean density 
in the region 11 - 17 arcminutes from the cluster center. 
In some instances, this was adjusted to avoid obvious 
over-densities in the background. 

In order to calculate R 30 , it was necessary to first de¬ 
termine the cluster ‘edge’, Rioo- This was difficult, as 
the edge was lost in the variations of the background 
density. In most clusters, a number of radii seemed ap¬ 
plicable as a possible Rioo and it was unclear which, if 
either, represents the ‘real’ cluster edge. 

The profiles in Figure O illustrate the difficulty in 
calculating R 30 particularly for the low Lx clusters. 
Aside from determining the cluster edge, the choice 
of background density can also have a large effect on 
R 30 . The unphysical nature of the very rapidly ris¬ 


Fig. 4.— The effect of varying the background density on the 
total number radius relations of MS 1455.0. The relations shown 
are calculated with the background density calculated as the mean 
density in the region 11 - 17 arcminutes from the cluster center 
with that density plus (stars) and minus (squares) its rms error. 

ing or falling total number-radius plots for some clus¬ 
ters ( e.g . RX0054.1) demonstrates this sensitivity. To 
further highlight this problem, Figure Q] shows how the 
total number of galaxies as a function of radius varies 
when the background density is changed. The back¬ 
ground level is changed by its ± 1 a error, resulting 
in wild variations in the total number profiles. This 
may be a particular problem for previous studies that 
have used a small field of view to image the clusters 
and then a separate global background estimate. The 
difficulty in determining R 30 Js further ill ustrated using 
t he c luster MS1455.0. ISmail. Edge. Ellis, fc Blandfordl 
(1998) calc ulated an R.m of 2. 2 arcmins for this clus¬ 
ter, whereas lFairlev etal . 12002} ) found 3.75 arcmins, and 
this study either 1.56 or 3.07 arcmins depending on the 
choice of the cluster edge. This could potentially be a 
significant problem for blue fraction studies, particularly 
when comparing samples from different authors, consid¬ 
ering that the blue fraction has been shown to depend 
on radius (}Butcher Jg_Oemler|1984tLAbrahain_eL-alji993 

Kodama fc Boweill200 ll lEUingson. Lin. Yee. fc Carlberel 

211010 - 

In an attempt to reduce some of the systematic un¬ 
certainties in calculating R 30 , it was decided to define 
another characteristic radius using the X-ray properties 
of this cluster sample. Assuming the cl uster core is viri- 
alised, in a simple self similar model ( Brvan fc Normanl 
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1 998; Eke. Navarro, fc FrenkSll9981 the virial radius, Ry, 
is given by 


Rv = 3.80/S?^A- ,/2 ( 1 + ^ ) - 3/a (ts|v) 1/2 '*» 1 - 

(1) 

with 

A, = (A c (lE,A)r> 0 )/(187T 2 fE) (2) 

taken from llArnaud. Aghanim. fc Neumann! (I2002f) . 
where A c (fl z ,A) is the density contrast, and (3t is the 
normalization of the vir ial relation GMy/2R.y = BrkT . 
A value of j3x = 1-05 (lEvrard. Metzler. fc Navarr 01199(11 
is used here. An analytic approxi mation for A c (fi z , A) is 
given in iBrvan fc Norman (tl998f) . 

Only three clusters in our sample have measured 
temperatur es to use in t he above rela tion: MS0353.6, 
M S1455.0 (IMushotzkv fc Scharfl 11997T) . and MS2137.3 
(lEttori Kr. Fa,hianlll999tl . Fbr the rest of the sample, it 
was necessary to convert the measured Lx to tempera¬ 
ture. This was done using the redshift dependent Lx—T 
relation, 


T = 


( L 


boi 


\B(l + z) A 


l/OL 


(3) 


where B = 3.11 ± 0.2 7 x 10 44 h 2 erg s 1 , A = 1.5 ± 
0.3, a = 2.64 ± 0.27 llMarkevitchlll994 iVikhlinin et alJ 
l 2 f) 02 ll - and I./,,,/ is the bolometric luminosity. 

A conversion from the ROSAT Lx to L boi was required 
before the above Lx — T relation could be used. This 
conversion was calculated using a model X-ray spectrum. 
The shape of the spectrum is dependent on the tempera¬ 
ture, and hence so is the conversion. The two steps were 
therefore combined iteratively, with an initial estimate 
of temperature used to calculate L boi , which was in turn 
used in the Lx — T relation to estimate the tempera¬ 
ture, which was again fed into the L boi conversion and so 
on until convergence. This procedure converged rapidly 
within 3-5 iterations. 

The complexity of the above method, combined with 
the existence of significant scatter in the Lx — T rela¬ 
tion, and the requirement of a virialised cluster core, call 
into question the reliability of Ry and its superiority over 
R 30 . Since all of these clusters show reasonably strong X- 
ray emission, the likelihood of the cores being virialised is 
high. To estimate how large an error both this conversion 
and the Lx — T scatter may have introd uced , the same 
method was applied to the data of IMushotzkv fc Scharll 
fl99l . They measured the temperatures of 38 clusters 
between a redshift 0.143 and 0.541 using ASCA. A typ¬ 
ical error of 20 % was found between the measured tem¬ 
peratures and those calculated with the Lx — T relation, 
which leads to only a 10% error in Ry. This error, and 
the previously mentioned uncertainties seem acceptable 
considering the particular advantages this method has in 
the removal of the fairly arbitrary decision on the po¬ 
sition of the cluster edge and its independence on the 
heavily fluctuating field density. The calculated values 
of L boi, T, and Ry are show in Table 0] 

Another advantage of X-ray data is revealed when de¬ 
termining the cluster center. The centroid of the X-ray 
emission yields the deepest point in the cluster poten¬ 
tial well. This often coincides with a bright elliptical 



Fig. 5. — B — R Color-Magnitude relations for the 5 lowest Lx 
clusters before (left) and after (right) K-corrections. The biweight 
fit to the k-corrected composite is shown as the solid line. 


galaxy. When this occurred, the center of this galaxy was 
used as the cluster center. Otherwise, the X-ray centroid 
was used. When only optical imaging data is present, 
the cluster center is often determined by the position of 
the br ightest cluste r galaxy. It has been demonstrated 
ijLvnam et al.112000 1 that this can often lead to an in¬ 
correct choice, potentially affecting the measurement of 
cluster properties. 

3.2. Statistical Field Correction 

In order to study the cluster galaxy population it was 
first necessary to correct for the presence of foreground 
and background field galaxies along the cluster line of 
sight. A field region was defined for each cluster as the 
area outside a radius of 3 Mpc. Having a large field of 
view allowed a field sample to be constructed for each 
cluster, which should better take into account variations 
in large-scale structure alon g the different lines of s ight. 

Following the method of iPimhhlet, et ~akl (I2002H the 
probability of each galaxy being a cluster member was 
determined in a color-magnitude grid of size 0.5 in mag¬ 
nitude and 0.1 in color by comparing the numbers of 
galaxies in the cluster and field regions. One hundred 
realisations of the cluster membership were then deter¬ 
mined using this probability grid. 

3.3. Fitting the Color-Magnitude Relation 

Although the color-magnitude distribution is domi¬ 
nated by the red galaxies defining the relation to be fit, 
there is still significant contamination from bluer galax¬ 
ies. To overcome this contamination, a robust biweight 
fit was made to the data, whic h assigns in c reasingly lower 
weights to outlying points. iBeers. Flvnn. fc Gebhardfl 
(19901 give scale (S bi) and location estimators (C bi) 
based on the biweight method. Using these estimators, 
the slope (M) and normalization (C) of the CMR are 
determined in turn. The CMR residual was defined as 
res = (M > 1 = magnitude + C) — color. First res was min¬ 
imized in Sbi with C = 0. This yielded the slope of 
the relation M. Then res was minimized in C bi with M 
fixed, giving the normalization of the relation. The pro¬ 
cedure was carried out for all clusters in each realization 
of the background subtraction, in both colors, within a 
radius of Ry/2, to a limiting magnitude My = -20. The 
final fit for each cluster was then determined by taking 
the median fit of the 100 realizations, and the error given 
as the rms of these fits. 

For the richer clusters, good fits were obtained to their 
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Fig. 6. — B — R Color-Magnitude relations for all 12 clusters. All galaxies to half a virial radius and to a limiting magnitude of My res i 
= -20 are shown. The large points are galaxies designated as cluster members and the small points as field members in one iteration of 
the statistical background subtraction. The solid near-horizontal line shows the median biweight fit to the relation and the shaded area 
the 3 sigma error on that fit from the background correction. For the five lowest Lx clusters the fit shown is that to the composite clus ter 
k-corrected back to the individual cluster redshift. The dashed vertical line shows the color dependent magnitude cut (see Section 13.61 . 
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TABLE 5 

B — R CMR FIT PARAMETERS AND COMPLETENESS. 
Cluster Slope Zero-Point Target Mag Comp 


M^ = -19 R- Band 


RX1633.6 -0.041 ± 

0.018 

3.17 ± 

0.35 

21.08 

23.6 

RX0333.0 -0.041 ± 

0.018 

3.19 ± 

0.35 

21.15 

23.7 

RX0210.4 -0.041 ± 

0.018 

3.28 ± 

0.35 

21.20 

23.5 

RX0054.0 -0.041 ± 

0.018 

3.33 ± 

0.35 

21.67 

23.4 

RX1606.7 -0.041 ± 

0.018 

3.37 ± 

0.35 

21.84 

23.4 

RX2237.0 -0.015 ± 

0.014 

2.82 ± 

0.26 

21.73 

23.3 

MS0407.2 -0.037 ± 

0.023 

3.04 ± 

0.43 

20.96 

22.9 

RX1418.5 -0.010 ± 

0.013 

2.76 ± 

0.26 

21.68 

23.3 

RX0256.5 -0.056 ± 

0.016 

3.74 ± 

0.32 

22.31 

23.7 

MS0353.6 -0.022 ± 

0.011 

2.95 ± 

0.22 

21.94 

23.3 

MS1455.0 -0.048 ± 

0.017 

3.34 ± 

0.33 

21.30 

22.7 

MS2137.3 -0.075 ± 

0.033 

4.01 ± 

0.67 

21.87 

23.3 


The median B — R CMR slope and zero-points calculated from the 
biweight fits to the 100 realizations of the background subtraction 
with the errors given as the rms of these fits and the R-band 
completeness limit estimated from the NMR. 

strong red sequences. However, in some of the poorer 
clusters, the sparsely populated red-sequences were not 
well defined and hence the fits were poor. To improve on 
this, a composite of the five lowest Lx clusters (Lx < 
10 44 ergs s _1 ) was made. The fit to this composite CMR 
would then be used for all of these clusters individually. 
To construct the composite, it was first necessary to k- 
correct the galaxy colors and magnitudes to account for 
the differences in redshift, and hence filter bandpass, be¬ 
tween the indi vidual clusters. This was d one using the K- 
corrections of iFukugita. Shimasaku. fe Ichikawal ( 1995 ?) 
with mean corrections of amplitude 0.08, 0.06 and 0.03 
in B, V and R respectively. Figure 0 shows 1 of the 
100 composite color-magnitude plots before and after K- 
corrections were applied, with the resultant biweight fit. 

Figure 0 shows B — R versus R color-magnitude dia¬ 
grams respectively for the 12 clusters. This filter combi¬ 
nation was chosen because it approximates U — V versus 
my at rest, thus spanning the 4000 A break and pro¬ 
viding a se nsitive measure of the age o f an old stella r 
population llDressler fc Shectmanlll987t iBruzual 1198.8 h 
The filled circles show galaxies selected as being in the 
cluster, and the dots show those in the field for this sta¬ 
tistical realization. In all plots, the median biweight fits 
of all the statistical realizations are shown as solid lines 
along with the 3 sigma errors on the fits as the shaded 
regions. Details of these fits are given in Table 0 For 
the low Lx clusters fitted as a composite, the zero-points 
have been k-corrected back to the redshift of the individ¬ 
ual clusters. We note that fits were also made to CMRs 
generated with the other filter combinations and that 
adding these to the analysis below does not change our 
results. 

By inspecting the color-magnitude plots, it can be seen 
that most of the clusters show clear CMRs. Those that 
do not, such as RX0333.0 and RX1633.6, are of low rich¬ 
ness as might be expected from their very low X-ray lu¬ 
minosities, and thus have very few galaxies with which 
to define a relation. When these clusters are combined 
together (Figure [ojl, a relation does become apparent. 
On plotting the fit from this composite cluster on the 
individual color-magnitude diagrams, a concentration of 
galaxies about the fit is then visible. 


3.4. The Dependence of the Color-Magnitude Relation 
on Cluster Mass 

As discussed in Section ^ the uniformity of color- 
magnitude relation provides tight constraints on any 
variation in the star formation history of early-type 
galaxies. This section will investigate whether or not 
the CMR is universal over the large range in cluster X- 
ray luminosity, and hence mass, present in the sample of 
clusters studied here. 

The properties of the CMR that will be considered here 
are the slope and zero-point, represented as the color at 
a fiducial magnitude. In order to investigate how these 
properties vary with cluster Lx , it is necessary to first 
correct for any dependence on redshift. Figured shows 
the slope and zero-point as a function of redshift for 
the clusters. The low Lx clusters, whose CMRs were 
fitted together in a composite, are shown as one point 
in these plots. Also shown on these plots are model 
predictions for the CMR using the PEG A SE2 stellar 
popu lation synthesis models llFioc fc R.occa-Volmerangdl 
flfifilt) n ormalised to the Com a CMR at a redshift 0.023 
(BowerJLucfi-V. J^ElliJ.iaa2bD . Clearly there is some dis¬ 
crepancy between the models and the fits is evident. 
Considering first the evolution of the slope, a large scat¬ 
ter is immediately noticeable, as well as systematically 
shallower slopes than the models predict. One possibil¬ 
ity for the steeper model slopes is the large apertures 
used to determine the colors. The Coma data, used to 
normalize the models, used a much smaller aperture to 
calculate the colors. Since there are radial color gradi¬ 
ents in elliptical galaxies, with redder colors nearer the 
centers, this has the effect of increasing the slope as only 
the inner regions are sampled in the brightest galaxies. 
Figure [7] also shows the zero-point of the CMR of each 
cluster, defined as the color of a galaxy lying on the fitted 
CMR that would have an absolute V band magnitude = 
-21.2 at a redshift zero. Again, there is a discrepancy be¬ 
tween the models and the fitted CMRs, with the cluster 
fits appearing redder than the models. The photometric 
checks discussed in Section ET51 suggest that the photom¬ 
etry is much more accurate than the 0.1 - 0.2 magnitude 
offset visible here. Allied to the consistent offset over 
both independent data sets, it appears that there may 
be a problem with t he model B magnitudes . This may 
not be too surprising. lEisenstein et all ( 2001 1 found that 
the PEGASE models were unable to reproduce the color 
evolution of Luminous Red Galaxies in the SDSS, par¬ 
ticularly in g — r, predicting either too red colors at low 
redshift or too blue at intermediate redshift. We have 
tied the models to the colors of the low z red galaxies in 
the Coma cluster which has made the models too blue at 
the intermediate redshift of the galaxies studied here. 

Although the model slopes are offset from the mea¬ 
sured slopes, the evolution with redshift appears to be 
consistent in form. This is seen in the lower left plots 
in Figure 0 where no gradient with redshift is apparent. 
Therefore, it is still possible to use the models to remove 
the redshift evolution of the slope in order to investigate 
any dependence on Lx , and this is plotted in Figure |H1 

For the zero-points, the models are clearly inadequate. 
Unlike the slopes, a tight relationship between zero-point 
and redshift is visible. This enables an empirical deter¬ 
mination of the relationship to be made and subtracted 
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Fig. 7.— The slope (left) and zero-point (right) of the B — R CMR as a function of redshift. The lines show the model predictions 
for different formation redshifts (see text for details). The bottom plots show the difference between the zp = 4 model and the measured 
slopes. The 5 lowest Lx clusters fitted as a composite are shown as one point (the dot inside a circle) at their mean redshift. 


in order to investigate the Lx dependence. It is safe to do 
this, as there is no dependence of Lx on redshift within 
this sample. 

Figure M shows the fits and residuals to the zero- 
point-redshift relations. The fits are simple, linear least- 
squares fits and provide good representations of the data. 
Figure ITHl shows the residuals to the fits as a function of 
X-ray luminosity. 

Inspecting Figures E and EH the color-magnitude re¬ 
lation shows no dependence on cluster X-ray luminosity 
either in its slope or zero-point over a factor of 100 in 
X-ray luminosity. This is confirmed by Spearman Rank 
tests, with correlation probabilities of 82% and 61% for 
the slope and zero-point respectively, consistent with no 
trend. The scatter in both the slope and zero-point Lx 
relations is consistent with being driven by the measure¬ 
ment errors. 

Even though the range in X-ray luminosity is large, 
the number of clusters considered is quite small, and the 
errors and scatter in these plots quite large. It is there¬ 
fore worth considering how much of a trend would still 
be consistent with these plots. To calculate this, a simple 
least-squares fit weighted by the errors was made to the 
relations in Figures E and EH A maximum variation of 
0.091 in the slope and 0.090 in the zero-point was allowed 
by the 2 sigma error to the fits over the entire Lx range 
of the clusters. The change in the slopes allowed is quite 
large, however, the maximum zero-point change is well 
constrained. 

If the change in zero-point is considered to be due to 
an age change, the maximum age change allowed can 
then be estimated. Using the PEGASE2 models, an age 
change of 1 Gyr corresponds to a color change of ~ 0.05 
in 13 — i? at z ~ 0.3. This means that an age difference of 
greater than ~ 2 Gyrs is ruled out at the 2 sigma level. It 
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Fig. 8.— The difference between the measured CMR and the 
model CMR slopes as a function of X-ray luminosity. Lx is given 
for the ROSAT band (0.4 - 2.0 keV) in units of 10 44 erg s -1 . 
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Fig. 9.— The CMR zero-point as a function of redshift (left) 
and the residual to the fits (right). 
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Fig. 10. — The difference between the measured and the model 
CMR zero-points as a function of X-ray luminosity (given for the 
ROSAT band (0.4 - 2.0 keV) in 10 44 ergs s -1 ). 

therefore appears that the global environment does not 
have a strong affect on the ages of early-type galaxies in 
clusters. 
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3.5. The Dependence of the Color-Magnitude Relation 
on Cluster Radius 

In the previous section we investigated how the CMR 
depends on the global environment of a galaxy cluster, as 
measured by its X-ray luminosity. Another environmen¬ 
tal change we can probe is the variation of local galaxy 
density as measured by cluster radius. One approach 
to investigating the CMR as a function of radius would 
be to try and fit the CMR in radial bins for each cluster. 
However, rapidly decreasing cluster galaxy densities, and 
thus increasing field subtraction errors, make this unfea¬ 
sible. Instead, all the clusters are combined together us¬ 
ing a m ethod similar to that described in lRimhblet et all 
(20021. Since we have already shown that the properties 
of red sequence galaxies are independent of the X-ray 
luminosity of the cluster in which they reside we can 
be confident that no significant biases due to the large 
range in L\ have been introduced when the entire clus¬ 
ter sample is combined together. The lowest Lx clusters, 
while contributing many fewer galaxies than the highest 
Lx clusters, do make a significant contribution, with the 
5 lowest Lx clusters together adding a similar number 
of cluster galaxies as one of the richest. Figure HD shows 
the surface density of the galaxies lying close to the CMR 
for each cluster. The clusters are generally centrally con¬ 
centrated with approximately circular density contours, 
and any irregularities should be smoothed out when the 
clusters are combined. 

Before combining, the CMR fits were subtracted from 
the galaxy colors of each cluster, i.e. for each galaxy the 
color difference (A {B — R)) from the fitted CMR was 
calculated. This was repeated for every galaxy in the 
cluster area (within 2 Ry) and for every galaxy in a 
combined field sample. The combined field sample was 
constructed by including all the galaxies in the cluster 
fields at a radius greater than 2.5 Ry from the cluster 
center. A combined sample is used here, rather then a 
field sample for each cluster as in the CMR fitting, both 
to allow a larger radius from the cluster center to be 


Mpc Mpc Mpc 

Fig. 11. — The galaxy surface density distributions. Only galax¬ 
ies with a color within 0.08 of both the fitted B — V and V — R color 
magnitude relations are considered. The positions are smoothed 
with a top hat filter of 500 Kpc in size. The lowest contour repre¬ 
sents 3 Mpc: “ 1 and each contour increases by 3 Mpc -1 . 

probed, and to maintain sufficient numbers to minimize 
the statistical errors in the field subtraction. 

For each cluster, the cluster and field samples were cut 
to a magnitude limit equivalent to a V absolute magni¬ 
tude of -20 at the cluster redshift, ~M* + 2. The cluster 
galaxies were then split into annuli of width one quarter 
of a virial radius, and binned into histograms of Acolor 
with bins of 0.1 magnitude. The field sample was simi¬ 
larly binned into a histogram. For each radial annulus, 
the field histogram was scaled according to the ratio of 
the areas and then subtracted from the cluster histogram. 
This is a much simpler method of field correction than 
that used when fitting the CMR, as the magnitude de¬ 
pendence has been removed, thus improving the number 
statistics in each bin, and enabling a more reliable back¬ 
ground subtraction over the same size color bin. The 
final step was to sum all the field subtracted histograms 
in each annulus to construct a composite of all the clus¬ 
ters. 

When combining a number of clusters, it is important 
to ensure the same sampling has been made for each clus¬ 
ter, both spatially and photometrically. The use of the 
virial radius ensured that the spatial scales of clusters 
were matched to each other. In order to compare the 
colors of galaxies at different redshifts, a K-correction is 
normally applied. This method was used when calcu¬ 
lating the ma gnitude limit for each cluster using the K - 
corrections of iFukugita. Shimasaku. fc Ichikawa! ( 1995 1. 
However, by calculating the color difference from the fit¬ 
ted CMR, the equivalent of a k and evolutionary correc¬ 
tion has been applied based on the empirical colors of the 
galaxies lying on the CMR of each cluster. This ‘correc¬ 
tion’ is likely to be more accurate than any determined 
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TABLE 6 

Gaussian fit parameters 
Radius Peak Sigma 

0 - 0.25.Rv 0.020 ± 0.009 0.159 ± 0.012 

0.25 R v - 0.5 R v -0.041 ± 0.012 0.142 ± 0.024 
0.5 Ry - 0.75 Ry -0.083 ± 0.035 0.184 ± 0.041 
0.75ijy - 1.25 R v -0.026 ± 0.133 0.139 ± 0.108 

by models, as it has been determined using the data on 
which it was to be applied. Its only limitation is that 
the equivalent of one correction was applied to all galax¬ 
ies regardless of how close their color lies to that of the 
CMR. However, as the galaxies of interest here are those 
lying close to the CMR, and the redshift range is small, 
any systematic errors would be much smaller than those 
introduced by using model K-corrections. 

Figure El shows the combined histograms, in 4 radial 
bins, for the 12 clusters. The solid histogram shows the 
field corrected data, and the dotted histogram shows all 
data prior to the field correction. Each panel shows the 
histogram in the radial bins of 0.25 Ry, except the final 
panel which shows the histogram in a larger bin, between 
0.75 and 1.25 Rv, in an attempt to increase the signal 
caused by the low cluster density at these radii. The 
central bin (top left panel) shows the clear peak of the 
cluster CMR. This gets smaller with increasing radii and 
disappears in the final bin (bottom right panel) at ~ Rv ■ 

To investigate any change in the color of the CMR with 
cluster radius, we wished to determine if and how the 
cluster peak color and width change. To do this, a Gaus¬ 
sian was fitted iteratively to the histogram data around 
the peak between -0.4 and 0.4 in A (B — R), thus ignoring 
the very bluest galaxies which become more significant at 
larger radii. The results of these fits are given in Table El 
and are plotted on Figure El as the dot-dashed curves, 
with the center of the Gaussian fit shown as a dashed 
line. To generate errors on the fits, random Gaussian 
noise scaled to each bin’s Poisson error was added to 
each histogram bin, and then the fit was recalculated. 
This was repeated 10000 times with the final fits and er¬ 
rors calculated as the mean and rms. A Gaussian fit was 
used to define the peak instead of other location estima¬ 
tors such as the biweight, as it could be fitted to data 
including negative bins, whereas the biweight cannot. 

The peak of the Gaussian fits show a clear trend of 
becoming bluer with increasing radius in the first 3 bins, 
with the fourth bin being very noisy with a large error 
on the fit. A total shift bluer of 0.10 ± 0.036 magni¬ 
tudes was observed over 0.75 Rv- The fact that the first 
bin shows a positive peak is reassuring, since the sub¬ 
tracted CMR was fitted within 0.5 Rv and so would lie 
between the first two radial bins considered here. The 
width of the CMR shows no evidence of any significant 
change with radius. However, the errors on the width are 
large, so definite conclusions regarding the CMR width 
are difficult. 

A clear variation in the color of the CMR with radius 
is seen for a combined sample of 12 intermediate redshift 
clusters. Removing the slope of the CMR before the ra¬ 
dial comparison is made ensures that any color change is 
independent of the magnitude distribution of the galaxies 
in a radial bin. This is important, as the magnitude dis¬ 
tribution will change, with the most massive, brightest 


and he nce reddest galax ies being located in the cluster 
center. iPimbblet, et, all ( 2002 1 measured the change in 
B — R. color as a function of radius in a combined sample 
of 11 rich X-ray selected clusters at a redshift of ~ 0.1. 
Since their clusters all had similar X-ray luminosities, 
they did not normalize the individual radii to a univer¬ 
sal scale, but just summed the clusters, measuring the 
radial color change in Mpc -1 . To directly compare the 
two trends, the virial radius of a cluster with the average 
X-ray luminosity and redshift of their sample was calcu¬ 
lated using the same method as for the clusters presented 
here, and the radii scaled accordingly. 

A comparison of the peak and width of the fitted 
Gaussians of the two studies is shown in Figure 1131 
Considering firstly the change in peak position, the 
two trends are similar, with the one presen te d her e 
being slightly steeper than IPimbblet et all ( 2002 1. 
Assuming that the color change refl ects a change in the 
age of the galaxies on the CMR. iPimbb let et alJl2002t 
Terlevich. Caldwell. fc Rowerl 120011 Mirahame^fl] 
1996), a steeper trend would be expected at a higher 
redshift due to the younger overall age of the galaxies, 
and thus the larger fractional age gradient with radius. 
Turning to the width of the fitted Gaussian, there 
is little evidence of any trend in the data p r esente d 
here, contrary to that found hv IPimbblet et aTl ( 2002 ). 

although the la rge err ors_make definite conclusions 

difficult (see also iBalogh et al.l (12004T )). 

Several other authors have reported similar radial 
color changes of the CMR. I Abraham et al.l ( 1996 1 stud¬ 
ied a single cluster, A2390, at a redshift of 0.23, 
with both photometry and spectroscopy, over a wide 
area. They found that the colors of galaxies on the 
CMR became bluer with radius with a gradient of 
d(g — r) / dlog io(r) = -0.079 out to 3.7 Mpc, equiva¬ 
lent to ~ 1.15 Rv calcula ted from the X-ray t emper ¬ 
ature meas ureme nt o f [Allen. *Et,tori. X Fa,bia.nl ( 200 11. 
iTerlevich. Caldwell & RoweiOToO 111 reported a bluen- 
ing of the U — V color magnitude relation in Coma as 
a function of cluster radius. Converting to the g — r 
color used by Abraham et al., they repo rt a change of 
d(q — r) / dloQin(r) — -0.024 ± 0.005. IPimbblet et all 
(20021 similarly converted their B — R color to g — r 
and find d{g — r)/dlogw{r) = -0.061 ± 0.011, which is 
intermediate between the other two studies. 

Converting the B — R color difference of the clus¬ 
ter sample presented here to g — r produces d(g — 
r) /dloq-\n(r) = -0.081 ± 0 .019, consistent with the trend 
of Abraham et all ( 1996 1 at the same redshift. Using 
simple evolutionary models presented in section 13.31 the 
age change corresponding to these color changes is calcu¬ 
lated, assuming that the galaxies are formed at a redshift 
> 2. This enables the rate of age change per log(r) to be 
determined. The results of these calculations are shown 
in Tabled 

All four studies produce very consistent results with 
a luminosity weighted age change of ~ 2.5 Gyrs per 
log(r), account i ng for the radial color cha n ge. Both 
lAbraham et all 1199611 and IPimbblet et al.l (1200211 ex¬ 
plain this change as a result of an increase in younger 
SO galaxies contam inating the CMR at larger radii. 
Illressler et a,ljfl997l l show that the fraction of SOs in 
clusters declines with increasing redshift suggesting that 
SO galaxies in clusters have been formed relatively re- 
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Fig. 12.— A (B — R) color histograms as a function of radius for the combined clusters. The solid histogram is field subtracted, and the 
dotted histogram includes the field galaxies. The dot-dash line shows the fitted Gaussian, and the dashed line shows the position of the 


Gaussian peak. 



TABLE 7 



Study 

z A(g-r)/log(r) 

A (g-r)/Gyr 

A(Age)/log(r) 

Terlevich et al. 

0.023 -0.024 ± 0.005 

0.0083 

-2.9 ± 0.6 

Pimbblet et al. 

0.12 -0.061 ± 0.011 

0.0151 

-2.5 ± 0.7 

Abraham et al. 

0.23 -0.079 ± 0.02 

0.0283 

-2.8 ± 0.7 

This Work 

0.28 -0.081 ± 0.019 

0.0323 

-2.5 ± 0.6 


Rate of change of g — r peak color as a function of log projected 
radius, the calculated change in color per Gyr of an old stellar 
population and the resulting rate of change of age as a function of 
projected radius from 4 photometric studies of clusters. 

cently. They further suggest that infalling spiral galaxies 
in the field are transformed into SO galaxies by the trun¬ 
cation of their star-formation. This may result in a gra¬ 
dient of younger SOs with cluster radius. iPimbhlet et all 
( 2002 1 suggest that the increasing width of the CMR his¬ 
togram they observe supports this, with a larger range 
of galaxy colors expected from the spread in age of the 
recently formed SOs. Unfortunately due to the large er¬ 
rors it is not possible to distin guish between n o trend and 
that of lPimbblet, et all (I2002T) here. Recently iHogg et all 


( 2004 1 have conducted a similar study of bulge domi¬ 
nated galaxies in the SDSS. They find evidence of a small 
shift in the peak color (< 0.02 mag g — r) and a broad¬ 
ening of the blue wing of galaxy color distribution with 
density further supporting this hypothesis. 

3.6. The Blue Fraction 

In the previous sections we have investigated the prop¬ 
erties of the red cluster galaxies as a function of cluster 
environment. We now turn to the blue cluster galaxy 
population and calculate th e fraction of blue galaxie s 
(ft) following the method of lButcher fc Oemled 11984T) . 
Butcher & Oemler defined the blue galaxies as those 0.2 
magnitudes bluer than the CMR in rest frame B — V. 
One approach to calculating ft, would therefore be to k- 
correct all of the data to the res t frame B—V. However, a 
simpler approach suggested bv lKodama fo Boweii ( 2001 1 
is to k-correct the blue galaxy criteria to the cluster red- 
shift. A B — V color 0.2 mag bluer than an elliptical 
galaxy falls between the colors of an Sab and Sbc at a 
redshift of zero. The blue fraction cut was thus defined 
as the color at the same position between a Sab and 
a Sbc at the cluster redshift, using the colors and K- 
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Fig. 13.— The CMR peak position (top) and width (bottom) 
as a fu nction of viria l radius for both this work (filled circles) and 
that of lPimbblet et aT H 200 I) (open circles). The peak position is 
shown relative to the central value. 

corrections of lFukugita. Shimasaku. fc Ichikawa! (1995). 
The blue fraction cuts are listed in Table 0 Butcher 
and Oemler also defined ff, to be measured to a limit¬ 
ing My of -20 at z = 0. At higher redshifts, this lim¬ 
iting magnitude depends on spec t ral ty pe. Again, the 
iFukugita. Shimasaku. fc Ichikawal 1 1 995 • colors and K- 
corrections were used to determine the limiting magni¬ 
tude for each cluster, which are shown as the dot-dot-dot- 
dash lines in Figure 0 These lines are close to vertical in 
most cases, due to the clusters being chosen at redshifts 
where the measurement filters map approximately onto 
rest UBV. The limiting magnitude for each cluster at the 
color of an elliptical galaxy is given in Table 0 

On calculating the blue fraction, a field correction was 
again required. However, since ft, is a simple measure 
which splits all the galaxies into just 2 groups, the com¬ 
plicated grid method was not necessary. Again, the back¬ 
ground was determined for each cluster at a radii > 3 
Mpc. The galaxies in both the cluster and held areas 
were then split into blue and red galaxies in the same 
manner. These numbers were then scaled by the rela¬ 
tive corrected areas, which are calculated in the same 
manner as in Section roi The held numbers were then 
subtracted from the cluster numbers to give fj, simply 
as residual blue galaxies divided by the residual total 
galaxy population. The error in ff, was determined using 
the Poisson errors on the numbers of galaxies in the blue 
and red, cluste r and held samples using the equations of 
iGehrelsi (119861) . 

When determining the held errors, it was important 
to consider the clustering of galaxies. To estimate this 
effect, the variance in the held region on the size of the 
cluster area was determined for each image. This was 
done by splitting the held into a grid, with each grid 


TABLE 8 
Blue fractions. 


Cluster Lim m# f& cut f& 


RX1633.6 

20.076 

0.182 

0.06 

RX0333.0 

20.145 

0.188 

-0.13 

RX0210.4 

20.201 

0.072 

0.21 

RX0054.0 

20.667 

0.093 

0.06 

RX1606.7 

20.839 

0.255 

-0.01 

RX2237.0 

20.729 

0.240 

0.06 

MS0407.2 

19.959 

0.172 

0.08 

RX1418.5 

20.678 

0.234 

0.10 

RX0256.5 

21.314 

0.303 

0.02 

MS0353.6 

20.937 

0.269 

0.04 

MS1455.0 

20.300 

0.200 

0.07 

MS2137.3 

20.869 

0.259 

0.02 


+0.39(0.29) 

- 0 . 29 ( 0 . 15 ) 

+ 0 . 63 ( 0 . 53 ) 

- 0 . 44 ( 0 . 26 ) 

+ 0 . 22 ( 0 . 20 ) 

- 0 . 16 ( 0 . 13 ) 

+ 0 . 19 ( 0 . 17 ) 

- 0 . 14 ( 0 . 09 ) 

+ 0 . 25 ( 0 . 20 ) 

- 0 . 20 ( 0 . 13 ) 

+ 0 . 24 ( 0 . 18 ) 

- 0 . 20 ( 0 . 13 ) 

+ 0 . 13 ( 0 . 11 ) 

- 0 . 11 ( 0 . 08 ) 

+ 0 . 11 ( 0 . 08 ) 

- 0 . 10 ( 0 . 06 ) 

+ 0 . 08 ( 0 . 06 ) 

- 0 . 07 ( 0 . 04 ) 

+ 0 . 09 ( 0 . 07 ) 

- 0 . 07 ( 0 . 05 ) 

+ 0 . 19 ( 0 . 12 ) 

- 0 . 17 ( 0 . 09 ) 

+ 0 . 13 ( 0 . 10 ) 

- 0 . 11 ( 0 . 07 ) 


Blue fractions for all clusters. The first error on each f 5 is that 
calculated using the rms background variation, and the second 
assuming a Poisson error in the background. 



Redshift 

Fig. 14.— Blue fraction as a function of redshift (left) and 
X-ray luminosity in units of 10 44 ergs s —1 (right). Error bars are 
calculated using the rms variation in the background. 


square having the same area as the cluster within which 
ff, is measured. The number of blue and red galaxies 
were calculated in the manner described above, and the 
mean and rms of this distribution determined. The rms 
was then used as the error in the held numbers as it 
represents how the Held beneath the cluster is typically 
varying, based on the local galaxy distribution. This 
rms error was always larger than the Poisson error on 
the whole background by a factor of ~ 8 , resulting in a 
increased error in ff, of ~ 50%. 

In order to investigate any dependence of the blue 
fraction on redshift and X-ray luminosity in a manner 
consistent with Butcher and Oemler, ff, was first calcu¬ 
lated within a radius of Ry /3 approximating R 30 . Table 
0 shows the resultant values of ff, in observed V — R 
(~ B — V at rest). Errors are given for purely Poisson 
variations on the number counts and for the rms back¬ 
ground variations. Figure m shows this data plotted as 
a function of both redshift and X-ray luminosity. It is 
clear that no obvious correlations exist and that there is 
a large scatter. Spearman Rank correlation tests confirm 
this. The large scatter and small range in redshift mean 
that it would be very difficult to draw any conclusions 
on the evolution of ft, with redshift from this data alone. 

This is not the case for the lack of a correlation of 
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Relative ram—pressure 

Fig. 15.— The blue fraction as a function of relative average 
ram-pressure. The ram-pressure is normalized to give the lowest 
Lx cluster RX1633.6 a value of 1. 

fb with Lx- This cluster sample covers a factor of ~ 
100 in X-ray luminosity with an even distribution of 
clusters throughout this range. The lack of any trend 
suggests that the blue fraction is independent of the 
global properties of the cluster, and that the distribu¬ 
tion of galaxy types is the same in all cluster cores re¬ 
gardless o f mass. This is in agreement with the work of 
iBalogh et al.l 6 2002 0. who show that the distributions of 
star formation rates in low Lx clusters at z ~ 0.25 are in¬ 
distinguishable from those of X-ray luminous clusters at 
similar redshifts, and lDe Pronris et al . (I2004IO who show 
that fb is independent of the velocity dispersion of local 
cl usters in the 2dF Ga laxy Redshift Survey. However, 
Margoniner. de Carvalho. Gal, & Diorgovskl (12001 1 and 
Goto et alJ (1200.3 ) both find a trend of increasing fj, with 
decreasing cluster richness. It should be noted that both 
studies use a fixed physical aperture within which they 
measure f&. This would correspond to a larger fraction of 
the virial radius, or R 30 , for the poor clusters, and could 
be driving this effect. 

It is worth briefly considering the implications of an 
independence of f& from Lx on the transformation mech¬ 
anisms discussed in Section El The most likely explana¬ 
tion for the Butcher-Oemler effect is that as time goes 
on there is a declining rate of accretion of field galax¬ 
ies on to clusters, combined with a declining star for¬ 
mation rate in these field galaxies, along with a mech¬ 
anism to tr uncate star formation w ithin the cluster en¬ 
vironment llKodama fc Bower 2001). If that truncation 
mechanism depends on the cluster mass then we might 
expect to see a correlation between ff, and Lx- For in¬ 
stance, if the truncation mechanism becomes ineffective 
for some mass threshold, all the accreted star forming 
galaxies would remain blue resulting in a higher blue frac¬ 
tion. For ram-pressure stripping we are able to directly 
determine the ram-pressure efficiency from the cluster 
Lx- Figure ESI shows the rest frame B — V blue frac¬ 
tion as a function of relative ram-pressure for all the 
clusters. The relative ram-pressure is calculated us¬ 
ing the relations P ram = picMV 2 al and V ga i oc kT$f 5 
and p 2 ICM oc Lx/kT^ L , where picm is the density of 
the intracluster medium, V ga i is the mean velocity of a 
galaxy through the cluster, and Tcl and Lx are intr¬ 
acluster medium temperature and X-ray luminosity re¬ 
spectively. The ram-pressure is normalized to give the 
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Fig. 16.— Blue fraction as a. function of Lx'- lo r 
4^Cu^_selectedcluster_samples_of IKodama & Bower (2001), 
|Smail J jE^g£ J _£llig J _^£lari^|Qr3 (1223) ■ aTlf l 

this work. The blue fractions from this work are calculated within 
Rv/3. 

lowest Lx cluster (RX0210.4) a value of 1. Figure ^3 
shows that the average ram-pressure varies by a factor 
of almost 40 in this sample, and yet no dependence with 
f h is observed. Further detailed nu merical simulations 
llStevens. Acreman. fc Ponmanl19991 ) and mo re sop histi- 
cated analytical models (|Fuiita fc Nagashimall999l) show 
that galaxies in the lowest Lx clusters observed here 
would not be expected to experience any ram-pressure 
stripping, even in the very central regions of the clus¬ 
ters where both their velocities and the cluster gas den¬ 
sity are highest. Therefore, it seems unlikely that ram- 
pressure stripping, or other mechanisms which depend 
on the global cluster properties such as tidal stripping 
or harassment, are the sole mechanism determining the 
blue fraction. 

Figure El shows the rest B — V bl ue fractions along 
with. 3 X-ra y selected c l uster sam ples (IKodama & Bowen 
2001|_ ISmail. E dge. El lis, fc Blandfordl Il998t 
Pim bble t et a . 120021 ) . This illustrates the exten¬ 
sion in the range of Lx sampled by this study, and 
demonstrates the lack of a correlation between f& and 
Lx over all these studies. 

A number of authors have reported that most 
cluste rs show increasing blu e fractions with ra¬ 
dius llB utcher fc Oemleil 11984 lAbraham et alJ Il996t 
iFairlev et all 120021 IKodama fc Bower! 120011) . This 
is consisten t with the morpholog y- densi ty relatio n 
fDressler. Oemler. Butcher, fc Gunnl 1994 k iDressleil 
11980) " and the observed increase i n the average star- 
formation rate with c l uster radius iTBijdp gh al] [11)99: 
Abraham et all Il996t iGomez et alJ 12003 iLewis et alJ 
2002 ). 

Figure E3 shows the average blue fraction calculated 
within a series of radii for the whole cluster sample. Val¬ 
ues of f;, within each radius are also calculated for 3 mag¬ 
nitude limits: the standard BO limit, and 1 magnitude 
either side of this. The trend of increasing ft with ra¬ 
dius is confirmed in this cluster sample with the main 
increase seeming to occur within Ry/2, with a near con¬ 
stant fraction beyond. The brightest galaxies show the 
largest trend, but the errors preclude making any quan¬ 
titative conclusion. This is not surprising, as the clus¬ 
ter core is expected to be dominated by bright elliptical 
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Fig. 17.— Blue fraction as a function of virial radius for the com¬ 
bined cluster sample. The blue fraction is calculated to three mag¬ 
nitude limits in My, -19 (crosses), -20 (circles) and -21 (squares) 

galaxies. 


Mv > -19 x 
Mv > -20 o 
Mv > -21 □ 



4. CONCLUSIONS 

We have analyzed the colors of bright galaxies in a 
sample of 12 X-ray selected clusters of galaxies at z ~ 
0.3. The clusters cover a large range in X-ray luminosity 
(10 43 - 10 45 erg s -1 ) and hence mass, with the lowest 
Lx clusters representing some of the least massive stud¬ 
ied at this redshift. Further more, the wide field of our 
imaging has allowed us to study the galaxy colors out 
to the virial radius of each cluster as determined by its 
X-ray luminosity, well beyond the core region typical of 
previous studies at this redshift. 

We have shown that the red galaxies in our clusters 
form a tight color-magnitude relation and measure its 
slope and zero-point. After removing the redshift evo¬ 
lution, we found that neither the slope nor zero-point 


depend strongly on cluster X-ray luminosity and hence 
mass. However, these constraints were limited by both 
the small number of clusters and the accuracy of the 
CMR fitting caused by the statistical field correction re¬ 
quired in purely photometric studied of galaxy clusters. 
By combining all the clusters together, we were able to 
measure the change in the zero-point of the CMR as a 
function of cluster radius. We found that the CMR be¬ 
came bluer by 0.1 ± 0.036 in B — R out to 0.75 II y- . If 
this was driven purely by an age change, it corresponds 
to a difference of ~ 2.5 Gyrs. 

We also measure the fraction of blue star forming 
galaxies as defined by Butcher and Oemler in our clus¬ 
ters. Again we find that there is no dependence of ff, on 
cluster X-ray luminosity and hence mass. However, we 
do again find a dependence on radius, particularly for 
the brightest galaxies. 

Overall, we can find no dependence of cluster galaxy 
population on the global cluster environment, i.e. clus¬ 
ter mass, but we do find a dependence on local environ¬ 
ment, i.e. galaxy density. The lack of any trends with 
cluster mass suggest that mechanisms like ram-pressure 
stripping, tidal stripping, or harassment, which depend 
strongly on cluster mass are unlikely to be solely respon¬ 
sible for changing the galaxy population from the blue 
star forming galaxies that dominate low density environ¬ 
ments to the red passive galaxies that dominate clusters. 
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